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Purpose: To investigate whether variants in a set of eight candidate genes are associated with diabetic retinopathy (DR) 
in a cohort of Chinese patients with type 2 diabetes mellitus (T2DM). 

Methods: Case-control study. Patients with T2DM were recruited from the Desheng community in urban Beijing and 
assigned into a DR group or diabetic without retinopathy (DWR) group, based on the duration of diabetes and grading 
of fundus images. Twenty-six single-nucleotide polymorphisms (SNPs) within eight candidate genes, including PPARy, 
vascular endothelial growth factor (VEGF) and its receptor kinase insert domain receptor (KDR), erythropoietin, aldose 
reductase, protein kinase C-[3, angiotensin-converting enzyme, and intercellular adhesion molecule 1, were analyzed 
using the MassARRAY genotyping system. 

Results: A total of 500 patients with T2DM (216 with DR and 284 with DWR) were enrolled in the study. Significant 
associations of DR were noted with genotypes of four SNPs— rs699947 (pO.OOl), rs833061 (p=0.001), rsl3207351 
(pO.OOl), and rs2146323 (p=0.006)— in the VEGF gene and one variant, rs2071559, in the KDR gene (p=0.034). 
After adjustment for covariates, significant association of DR remained with the homozygous genotype of the minor 
allele for the SNPs rs699947 (odds ratio [OR] = 3.54, 95% confidence interval [CI]: 1.12-11.19), rs833061 (OR = 3.72, 
95% CI: 1.17-11.85), rsl3207351 (OR = 3.76, 95% CI: 1.21-11.71), and rs2146323 (OR = 2.8, 95% CI: 1.46-5.37) in the 
VEGF gene as well as the SNP rs2071559 (OR = 1.62, 95% CI: 1.08-2.41) in the KDR gene. However, only rs699947 
and rsl3207351 in the VEGF gene remained statistically significant after Bonferroni correction. No associations were 
found in other genes tested. 

Conclusions: These data expanded previous observations on the association of DR with variants in the VEGF gene in 
Chinese patients with T2DM. Moreover, a possible association between DR and KDR polymorphisms is suggested. 



Diabetic retinopathy (DR) is one of the most common 
complications of diabetes and the leading cause of irrevers- 
ible visual loss in people of working age in the United States 
[1]. The etiology of DR is complex and largely unknown. 
Hyperglycemia and long duration of diabetes along with 
hypertension are widely recognized as major risk factors for 
its development [2,3]; however, these factors seem to explain 
only a part of the risk of DR [4,5]. In clinical practice, for 
example, some patients with poorly controlled diabetes or a 
long duration of diabetes do not develop retinopathy, whereas 
others with relatively good glycemic control have advanced 
retinopathy [6]. Accumulated evidence has suggested that 
genetic factors may play a role in the development and 
progression of DR, with heritability estimated to be as high as 
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27% [7] for DR and 50% for proliferative diabetic retinopathy 
(PDR) [8]. 

Several studies have used a genome-wide association 
approach [9-12] and multipoint sib-pair analysis [13,14] to 
detect regions of the genome that potentially contain genes 
involved in the etiology of DR. While different regions of the 
genome have been implicated, the reported regions were not 
consistently identified among the studies. In addition, many 
studies have taken a candidate gene approach to investigate 
the genetic etiology of DR, implicating the potential candidate 
genes in DR etiology, which include genes for PPARy [15], 
vascular endothelial growth factor (VEGF) [16-19], erythro- 
poietin (EPO) [20,21], aldose reductase (AKR1B1) [22-25], 
protein kinase C (PKC)-|3 [26,27], angiotensin-converting 
enzyme (ACE) [28], and intercellular adhesion molecule 1 
(ICAM-1) [29-31]. The validity of these results, however, 
remains unconfirmed, with few of the reported susceptibility 
variants having been consistently replicated [25]. 
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Kinase insert domain receptor (KDR), also called VEGF 
receptor (VEGFR) 2, is believed to be responsible for the 
majority of the angiogenic and permeability-enhancing 
effects of VEGF [32,33]. An increase expression of KDR has 
been reported in the retinas from diabetic rats [34,35]. One 
study investigated the fibrovascular tissues obtained from 
22 PDR patients and suggested the co-expression of KDR 
and neuropilin-1 facilitated fibrovascular proliferation in 
DR [36]. Moreover, microvascular expression of KDR was 
found to be associated with leaky vessels in diabetic retina 
from eyes of diabetic donors as compared to the retina from 
eyes of nondiabetic control donors [37]. Taken together, these 
findings imply that KDR might be a major contributor in the 
pathogenesis of DR. Furthermore, the KDR gene has been 
found to be associated with other pathogenic conditions 
featured for vascular permeability and angiogenesis, such as 
age-related macular degeneration (AMD) [38], cancer [39], 
and coronary heart/artery disease [40,41]. However, genetic 
association studies with DR have not been performed for the 
KDR gene. In this study, we tested whether the candidate 
genes, including PPARy, VEGF, EPO, AKR1B1, PKC-|3, 
ACE, and ICAM-1 genes as well as the KDR gene, were 
associated with DR in our independent cohort of Chinese 
patients with T2DM. 

METHODS 

Subjects and clinical evaluation: The study protocol was 
approved by the Ethics Committee of Beijing Tongren 
Hospital (TRECKY200907) and adhered to the tenets of 
the Declaration of Helsinki. Written informed consent 
was obtained from all participants before their enrollment. 
Patients with T2DM were recruited between November 2009 
and September 2011 from the Desheng Community of urban 
Beijing. Diabetes was defined as a self-reported history of 
physician-diagnosed T2DM treated with insulin, oral hypo- 
glycemic agents, or diet only; or by a fasting plasma glucose 
(FPG) concentration of 7.0 mmol/1 (126 mg/dl) or more in 
at least two previous examinations; or a random plasma 
glucose concentration of >11.1 mmol/1 (200 mg/dl). All 
subjects underwent a standardized evaluation consisting of a 
questionnaire, ocular and anthropometric examinations, and 
a laboratory test. The questionnaire elicited basic information 
(age, sex, ethnicity, income, education), lifestyle information 
(such as smoking and alcohol intake), health status informa- 
tion (such as the use of insulin therapy and any history of 
systemic diseases), and family history of diseases. Anthro- 
pometric parameters included body weight and height, waist 
and hip circumference, and three measurements, 5 min apart, 
of blood pressure in a resting state. Body mass index (BMI, 
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kg/m 2 ) was calculated according to the height and weight 
of the participant, and waist-to-hip ratio was calculated. 
A comprehensive ophthalmological examination included 
corrected visual acuity, slit-lamp biomicroscopy, and dilated 
fundus photography. Seven fields of 30°color fundus photo- 
graphs with stereoscopic images of the optic disc and macula 
were taken through the dilated pupils of each patient, using 
a digital fundus camera (Zeiss Visucam Pro; Oberkochen, 
Germany). 

Based on the duration of diabetes and grading of fundus 
photographs, patients were assigned to the diabetic-without 
retinopathy (DWR) group if they had more than 10 years of 
T2DM with no signs of DR (microaneurysms, hemorrhages, 
exudates) or if they had more than 15 years of T2DM with 
fewer than five microaneurysms. Patients with five or more 
microaneurysms in at least one eye were assigned to the DR 
group. Patients who did not meet the DWR or DR criteria 
were excluded from this study. The duration of diabetes was 
defined as the interval between the first diagnosis of diabetes 
and the time of enrollment in the present study. 

Laboratory assays: Fasting venous blood samples were 
collected for measurement of FPG, glycosylated hemoglobin 
Ale (HbAlc), creatinine, uric acid, and lipid profile (levels 
of total cholesterol, triglycerides, and high-density and low- 
density lipoprotein cholesterol), which were measured in 
an automated system with reagents for routine biomarkers. 
HbAlc was assessed by the enzymatic method using a 
Hitachi analyzer 7080 (Hitachi, Ibaraki, Japan). A first-void, 
midstream, morning, spot urine sample was collected, and 
albuminuria was measured by immunonephelometry with a 
Roche/Cobas C501 analyzer (Ibaraki). High albuminuria was 
defined as >20 mg/1. 

Candidate genes and single-nucleotide polymorphisms selec- 
tion: Eight candidate genes that are involved in the patho- 
genesis of DR or previously identified as being associated 
with DR were tested, including the PPARy, VEGF and its 
receptor (KDR), EPO, AKR1B1, PKC-P, ACE and ICAM-1 
genes. Twenty-six single-nucleotide polymorphisms (SNPs) 
in these genes were chosen initially for availability of assays, 
informativeness, and spacing across the gene. More informa- 
tion about the candidate genes is given in Table 1. It should 
be noted that the VEGF gene has been reported in previous 
literature [16], and the present study repeated the associa- 
tion between the VEGF gene and DR in an expanded larger 
sample size. 

Genotyping: Blood samples were collected from all partici- 
pants and stored at -80 °C before DNA extraction. Genomic 
DNA was extracted from venous blood leukocytes using a 
genomic DNA extraction and purification kit (TIANamp 
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Table 1. Description of candidate genes investigated 



Gene 



Gene function 



Locus 



References 



PPARy 
KDR 

VEGF 

EPO 
AKR1B1 

PKC-P 

ACE 

ICAM-1 



responsible for differentiation of fibroblasts to adipo- 
cytes and regulation of their function 

functions as the main mediator of VEGF-induced 
endothelial permeability, proliferation, survival, 
migration, tubular morphogenesis and sprouting 

a critical angiogenic and vasopermeability factor that 
is upregulated in various conditions involving retinal 
ischemia 

been linked to angiogenesis, vasculogenesis, and 
neuroprotection 

the first and rate-limiting enzyme in the polyol 
pathway 

PKC family members phosphorylate a wide variety 
of protein targets and are known to be involved in 
diverse cellular signaling pathways 

involved in the conversion of angiotensin I to angio- 
tensin II (ATII). 

a member of the immunoglobulin superfamily of 
adhesion molecules, involved in blood-retina barrier 
breakdown, capillary nonperfusion and endothelial 
cell damage and death 



3p25 
4qll-ql2 

6pl2 

7q21-7q22 
7q35 

16pll.2 

17q23 

19pl3 



Maciej et al. 2008 [13] 



Galan et al. 2010 [46] 



Churchill et al. 2008 [15] Al-Kateb et 
al. 2007 [16] 

Tong et al. 2008 [18] Abhary et al. 2010 
[19] 

Richeti, F et al. 2007 [22] Abhary et al. 
2009 [23] 

Uthra et al. 2010 [24] Ikeda et al. 2004 
[25] 

Kondo et al. 2009 [26] 



Kamiuchi, et al. 2002 [21] L. Liu, et al. 
2006 [22] Petrovic, et al. 2008 [23] 



Locus: the chromosomal region of gene; p: long arm of chromosome; q: short arm of chromosome 



Swab DNA Kit; Tiangen Biotech, Beijing, China). Study 
participants were genotyped for the SNPs using Sequenom 
MassARRAY technology (Bioyong Technologies, Beijing, 
China). All DNA samples passing initial quality checks were 
plated at a concentration of >5 ng/ul for processing on the 
platform. Quality measures taken into account for genotyped 
SNPs to be excluded from the subsequent analysis were 
minor allele frequency (MAF) <0.05, genotyping success 
<80%, and failed Hardy-Weinberg equilibrium (HWE) test 
in control samples (p<0.001). Quality control criteria for each 
SNP were implemented in our data set. 

Statistical analysis: Statistical analysis was performed using 
the R statistical analysis package. Baseline characteristics of 
diabetic patients in the DR and DWR groups were compared 
using a t test for continuous variables or a chi-square test 
for categorical variables. The genotype frequencies of genes 
were checked for HWE in all groups using a chi-square test. 
The chi-square test was also used to analyze the distribution 
of genotypes and alleles. When the expected frequency was 
less than 5, Fisher's exact test was used. The pairwise linkage 
disequilibrium was calculated using Haploview, version 4.2 
(Broad Institute, Cambridge, MA). Univariate and multivar- 
iate logistic regression analyses were performed to identify 
the association between SNPs/haplotypes and the presence of 
DR. Statistical results were expressed as p values, odds ratios 



(OR), and 95% confidence intervals (CI). The Bonferroni 
correction (p<0.05 divided by the number of SNPs analyzed) 
was used to account for multiple testing. 

RESULTS 

In total, 1,433 subjects with T2DM were enrolled from 
the community, of which, 500 subjects (284 in the DWR 
group and 216 in the DR group) participated in the genetic 
study based on the inclusion criteria. The distribution of 
demographic, behavioral, and clinical characteristics of 
participants are listed in Table 2. In the DR group, 59 patients 
were further classified as PDR. Compared with the DWR 
group, individuals in the DR group had a younger age of 
diabetic onset (p=0.008); higher percentage of high albumin- 
uria (p<0.001); higher levels of creatinine (p=0.027), FPG 
(p<0.001), and HbAlc (p<0.001); and were more likely to use 
insulin (p<0.001). Duration of diabetes in the DWR group 
was longer than the DR group (p=0.013), which could be due 
to the definition of DWR in this study (more than 10 years 
of T2DM without retinopathy). No statistically significant 
differences were found between the DR and DWR groups 
in terms of gender, BMI, waist-to-hip ratio, blood pressure, 
serum levels of uric acid, or lipid profile. 
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Table 2. Clinical and biochemical markers for the studied croups 



Clinical characteristics 


DWR (n=284) 


DR (n=216) 


p value 


Age of diabetic onset (years) 


52.7±7.61 


50.67±9.41 


0.008 


Sex (Male/Female) 


111/173 


103/113 


0.054 


Duration of diabetes (years) 


14.79±4.97 


13.44±7.22 


0.013 


BMI (kg/m 2 ) 


25.26±3.95 


25.8±4.13 


0.14 


WHR 


0.92±0.06 


0.93±0.06 


0.5 


High albuminuria(-/+) 


240/38 


143/68 


<0.001 


Systolic blood pressure (mmHg) 


136.7±16.65 


138.5±17.18 


0.22 


Diastolic blood pressure (mmHg) 


77.65±9.36 


79.23±9.69 


0.066 


Insulin therapy (yes/no) 


92/191 


118/97 


<0.001 


HbAlc (%) 


6.97±1.40 


7.84±1.70 


<0.001 


FPG (mmol/1) 


8.09±2.43 


9.22±3.21 


<0.001 


Creatinine (umol/1) 


67.43±17.69 


73.66±49.84 


0.027 


Uric acid (umol/1) 


281.9±82.42 


281.1±76.58 


0.98 


Cholesterol (mmol/1) 


5.05±0.98 


5.18±1.17 


0.19 


Triglycerides (mmol/1) 


1.55±0.93 


1.73±1.52 


0.11 


HDL cholesterol (mmol/1) 


1.25±0.30 


1.23±0.30 


0.38 


LDL cholesterol (mmol/1) 


3.07±0.83 


3.08±0.90 


0.83 



The data are expressed as mean ± standard deviation (SD) or number of subjects. Differences between diabetic retinopathy (DR) and 
diabetic without retinopathy (DWR) groups were compared using / test or Chi-square test. BMI, body mass index; WHR, waist and hip 
ratio; HbAlc, glycosylated hemoglobin Ale; FPG, fasting plasma glucose; HDL, high-density lipoprotein; LDL, low-density lipoprotein. 



Basic information of SNPs analyzed with allelic distri- 
butions is shown in Table 3. Out of the 26 SNPs selected, 
rsl 1645239 in the PKC-P gene had a low genotyping call rate 
(<52%); rsl805192 in the PPARy gene and rsl799969 in the 
ICAM-1 gene had a minor allele frequency of less than 5% in 
the studied population. These three SNPS were not studied 
further. The remaining 23 SNPs were tested in the DR and 
DWR groups for any departure from HWE. All were in HWE 
except rs3900008 (p<0.001), which failed to be genotyped in 
53 individuals, indicating a possible genotyping error. The 
SNP rs3900008 was therefore excluded, and a total of 22 
SNPs were included for subsequent analysis. For each of the 
SNPs evaluated, the number of participants varied slightly 
depending on the genotype call rate (>94%). 

Table 4 shows the genotypic distribution in the DR 
and DWR groups for the 22 informative SNPs. A strong 
association with DR was detected at three promoter vari- 
ants (rs833061, rsl3207351, rs699947) and one intronic SNP 
(rs2146323) in the VEGF gene, with the minor alleles of SNPs 
rs833061, rsl3207351, rs699947, and rs2146323 showing a 
similar magnitude of risk (OR>1.52; p<0.003). The ORs for 
the homozygote minor alleles were 3.80 (95% CI: 1.91-7.57) 
for rs699947, 3.6 (95% CI: 1.79- 7.22) for rs833061, 3.7 (95% 
CI: 1.86-7.39) for rsl3207351, and 2.8 (95% CI: 1.46-5.37) 



for the intronic SNP rs2146323. The three promoter variants 
rs833061, rsl3207351, rs699947 in the VEGF gene showed 
strong linkage disequilibrium (r 2 >0.98). SNP rs2071559 in the 
KDR gene was also significantly associated with DR with an 
OR of 1.41 (95% CI: 1.07-1.86) for the minor allele and 1.57 
(95% CI: 1.05-2.28) for the minor homozygote. 

Table 5 presents the multivariate analyses data adjusted 
for variables including the age of diabetic onset, duration of 
diabetes, high albuminuria, insulin use, HbAlc, and creati- 
nine levels. The association with DR remained statistically 
significant with the homozygous genotype of the minor allele 
for SNPs rs699947 (OR=3.74, 95% CI: 1.75-8.00), rs833061 
(OR=3.50, 95% CI: 1.62-7.54), rsl3207351 (OR=3.64, 95% 
CI: 1.71-7.77), and rs2146323 (OR=2.88, 95% CI: 1.4-5.93) in 
the VEGF gene as well as SNP rs2071559 (OR=1.88, 95% 
CI: 1.01-3.48) in the KDR gene. However, only rs699947 and 
rsl3207351 in the VEGF gene remained statistically signifi- 
cant after Bonferroni correction for the 22 SNPs analyzed. 
No statistically significant association with DR was found 
for all other SNPs in genes for EPO, AKR1B1, PKC-p 1 , ACE, 
and ICAM-1. 

Haplotype analysis revealed five SNP blocks (Figure 
1). Only two haplotypes in the VEGF gene showed strong 
evidence of an association with DR, and the others did not 
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Table 3. Basic information of Alleles and allele association results Between DR and DWR croups 



Gene 



locus 



SNP 


position 


Minor 


MAF 


HWpval 


OR (CI) 


P 


Allele 


DR DWR 


value 


rsl805192 


12,396,238 


G 


0 










rs2071559 


55,687,123 


C 


0.34 


0.39 


0.022 


1.41(1.07,1.86) 


0.011 


rsl3109660 


55,665,437 


A 


0.31 


0.63 


0.49 


1.06(0.8,1.41) 


0.66 


rsl870378 


55,661,210 


A 


0.48 


0.67 


0.33 


0.92(0.71,1.21) 


0.55 


rsl870377 


55,667,731 


A 


0.48 


0.4 


0.73 


0.95(0.73,1.23) 


0.67 


rs699947 


43,844,367 


A 


0.28 


0.009 


0.51 


1.63(1.22,2.18) 


0.001 


rs833061 


43,845,464 


C 


0.28 


0.021 


0.33 


1.53(1.14,2.04) 


0.003 


rsl3207351 


43,845,772 


A 


0.28 


0.01 


0.41 


1.58(1.18,2.12) 


0.001 


rs2010963 


43,846,328 


C 


0.44 


0.41 


0.9 


0.87(0.67,1.13) 


0.29 


rs833069 


43,850,557 


G 


0.44 


1 


0.9 


0.92(0.71,1.2) 


0.53 


rs2146323 


43,853,073 


A 


0.26 


0.003 


0.49 


1.52(1.13,2.05) 


0.004 


rs3025021 


43,857,141 


T 


0.16 


0.58 


1 


1.36(0.95,1.97) 


0.083 


rs3025039 


43,860,514 


T 


0.09 


1 


0.38 


1.34(0.96,1.88) 


0.077 


rs507392 


100,157,872 


C 


0.2 


0.4 


0.03 


0.93(0.67,1.28) 


0.64 


rs551238 


100,159,464 


C 


0.2 


0.39 


0.045 


1.04(0.75,1.45) 


0.8 


rsl617640 


100,155,234 


G 


0.19 


0.1 


0.14 


1.18(0.85,1.66) 


0.31 


rs759853 


133,794,498 


A 


0.07 


0.79 


0.34 


1.31(0.92,1.87) 


0.11 


rs3900007 


23,754,453 


G 


0.42 


0.49 


0.22 


1.01(0.78,1.31) 


0.95 


rs3900008 


23,754,514 


C 


0.3 


<0.001 


<0.001 


1.05(0.78,1.41) 


0.75 


rs 11645239 


23,754,563 


G 


0.36 










rs2575390 


23,754,255 


G 


0.06 


1 


1 


0.95(0.53,1.69) 


0.85 


rs3760106 


23,753,297 


A 


0.06 


1 


1 


1.02(0.57,1.84) 


0.95 


rsl800764 


58,904,261 


C 


0.39 


1 


1 


1.23(0.94,1.6) 


0.12 


rs9896208 


58,929,841 


T 


0.42 


0.26 


0.021 


0.95(0.73,1,24) 


0.71 


rs5498 


10,256,683 


G 


0.31 


0.64 


0.1 


1.32(1,1.75) 


0.044 


rsl799969 


10,255,792 


A 


0.003 











PPARy 



KDR 



VEGF 



EPO 



AKR1B1 



PKC-P 



ACE 



ICAM-1 



16 



17 



19 



MAF: minor allele frequency; HWpval: p value for Hardy-Weinberg (H-W) equilibrium; OR: odds ratio; CI:95% confidence interval; 
Genotype distributions for SNPs were in Hardy-Weinberg (H-W) equilibrium, except rs3900008 in both groups. 



show any statistically significant association (Table 6). The 
frequency of haplotype ACGA in the VEGF gene, defined 
by the minor alleles of the three SNPs rs699947, rs833061, 
rsl3207351 at the promoter region and rs2010963 at the 5' 
untranslated region, was 32.4% in the DR group and 23.2% 
in the DWR group (p=0.002, OR=1.58, 95% CI: 1.19-2.10). 
The haplotype AA defined by the minor alleles of the SNPs 
rs833069 and rs2146323 at intron 2 of the VEGF gene 
was 30.7% in the DR group and 22.5% in the DWR group 
(p=0.004), with an OR of 1.53 (95% CI: 1.15-2.04) for the 
increased risk of DR. 



DISCUSSION 

We analyzed the possible association between polymorphisms 
in eight candidate genes and DR in a well-defined cohort of 
Chinese patients with T2DM. Our data showed significant 
association between SNPs in the VEGF gene and the risk of 
DR after adjustment for possible confounding risk factors. 
However, no statistically significant association was found for 
other genes, including KDR, PPARy, EPO, AKR1B1, PKC-P, 
ACE, and ICAM-1. 

The pathogenesis of DR is complex with multifactorial 
biochemical causes influenced by genetic and environmental 
factors. Many vasoactive factors stimulated by hypergly- 
cemia or oxidative stress have been identified [1]. Among 
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Figure 1. Linkage disequilibrium plot generated by Haploview 4.2 software. Five haplotype blocks (bold) were identified for single-nucleotide 
polymorphisms (SNPs) in the kinase insert domain receptor, vascular endothelial growth factor gene, protein kinase C-(3 gene, and in eryth- 
ropoietin gene. Linkage disequilibrium (LD) is displayed as the pairwise D' value multiplied by 100 and given for each SNP combination. 
Shading represents the magnitude and significance of the pairwise LD, with a red-to-white gradient reflecting higher-to-lower LD values. 
Red diamond without a number corresponds to a D' value of 1.0. 



them, VEGF has been characterized as a critical angiogenic 
and vasopermeability factor that is upregulated in various 
retinal ischemic disorders, such as retinopathy of prematurity, 
AMD, and DR [32], whereas anti-VEGF therapy has been 
shown to be effective for patients with AMD [42] or diabetic 
macular edema [43]. The human VEGF gene is highly poly- 
morphic and has a large number of individual SNPs examined 
in relationship with DR [16-19,44]. In our previous study, 
three SNPs (rs833061, rsl3207351, rs699947) at the promoter 
region of the VEGF gene showed a strong association with 
the risk of DR [16]. This current study with an expanded 
sample size confirmed our previous report. Discussions of 
the three promoter SNPs have been described in detail in 
our previous literature [16]. Moreover, the association of DR 
with an intronic SNP rs2 146323 in the VEGF gene has been 
suggested in this current study. Although the intronic region 
of VEGF is believed to contain binding sites for transcrip- 
tion factors that regulate VEGF production [45], there is no 
functional data to date on the intronic SNPs. 

There are three high-affinity VEGF tyrosine kinase 
receptors: VEGFR-1, VEGFR-2 (KDR), and VEGFR-3 that 
are expressed mainly in the vascular endothelial cells [32,33]. 
KDR is the primary responder to the VEGF signal regulating 
endothelial cell migration and proliferation [32,33]. Several 
lines of evidence have implicated that KDR plays a role in the 
etiology of DR [34-37]. Although it became insignificant after 
Bonferroni correction, the present study suggested a possible 
association between the SNP rs2071559 in the KDR gene and 
DR. Further studies with a larger sample size to confirm this 



observation would be warranted. The present study failed 
to find an association of DR with other previously reported 
genes, including ACE, AKR1B1, PKC-|3, PPARy, ICAM-1, 
and EPO. Discrepancies between our study and previous 
reports could be due to different ethnic origin, diagnostic 
criteria, or sample size. 

Strengths of the current study include the enrollment of 
participants from a relatively homogeneous population of a 
single community, collection of detailed information on a 
variety of risk factors, and the standardized fundus photo- 
graph and assessment of DR severity. Additionally, previous 
genetic association studies for DR have been limited to studies 
of one or a modest number of candidate genes [25,46]. To our 
knowledge, this is the first and largest study in the Chinese 
population to address the question of whether multiple candi- 
date genetic polymorphisms are associated with DR. 

In summary, the present study confirmed our previous 
observation on the association between DR and VEGF 
polymorphisms independent of other clinical factors in the 
Chinese population. Moreover, a possible association between 
DR and KDR polymorphisms is suggested and further studies 
are warranted. 
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Table 6. Haplotype association analysis between polymorphisms and risk of 
diabetic retinopathy in patients with type 2 diabetes mellitus 



Haplotype blocks 


DR group (%) 


DWR group (%) 


OR (95%CI) 


Chi Square 


p value 


Block 1 












AAG 


48.7 


47.2 


1.07 (0.82, 1.39) 


0.26 


61 


TGA 


29.5 


31.2 


0.93 (0.69, 1.23) 


0.28 


0.59 


TGG 


20.7 


20.8 


0.99(0.72-1.37) 


0.001 


0.98 


Block 2 












ACGA 


32.4 


23.2 


1.58 (1.19,2.10) 


9.87 


0.002 


CTCG 


42.1 


44.6 


0.90(0.69, 1.16) 


0.68 


0.41 


CTGG 


25.2 


31.6 


0.73(0.55, 0.97) 


4.8 


0.03 


Block 3 












AA 


30.7 


22.5 


1.53 (1.15,2.04) 


8.39 


0.04 


AC 


25.8 


32.4 


0.72 (0.55, 0.96) 


5.16 


0.02 


GC 


43.4 


44.8 


0.94(0.73, 1.21) 


0.23 


0.63 


Block 4 












CGA 


57.3 


57.9 


0.97(0.75, 1.26) 


0.044 


0.83 


CGG 


37.1 


36.3 


1.04(0.80, 1.35) 


0.068 


0.79 


GAG 


5.6 


5.8 


0.97(0.57, 1.68) 


0.009 


0.92 


Block 5 












CCG 


17.5 


20 


0.86(0.62, 1.20) 


0.79 


0.37 


TAT 


80.6 


79.4 


1.16(0.84, 1.61) 


0.79 


0.37 



Data were given as frequency of each haplotype within diabetic retinopathy (DR) group or diabetic without retinopathy (DWR) group. 
Block 1 included rsl870378, rsl3109660, and rsl870377 in kinase insert domain receptor gene, block 2 included promoter SNPs 
rs699947, rs833061, rsl3207351 and rs2010963 at 5' UTR, and Block3 included SNPs two intronic SNPs of rs833069 rs2146323 in 
vascular endothelial growth factor gene, Block4 included SNPs rs3760106 rs2575390 rs3900007 in protein kinase C-(3 gene, block 5, 
rsl617640 rs507392 rs551238 in erythropoietin gene. Only haplotypes with frequency over 3% were included in this table. OR indicates 
odds ratio and CI refers to confidence interval. 
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